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Abstract 

We perform, with the help of cloud computing resources, extensive Langevin simulations which 



i-G ■ provide compelling evidence in favor of a general markovian framework for unbiased polymer 

■ ' translocation. Our statistical analysis consists of careful evaluations of (i) two-point correlation 

functions of the translocation coordinate and (ii) the empirical probabilities of complete polymer 
a ' 
c/3 ■ translocation (taken as a function of the initial number of monomers on a given side of the mem- 

C^ . brane) . We find good agreement with predictions derived from the Markov chain approach recently 

addressed in the literature by the present authors. 
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I. INTRODUCTION 



The phenomenon of polymer translocation through membrane pores has been the subject 
of a fiury of research activity in recent years [l|, a fact related to its relevance to the under- 
standing and development of important biotechnological processes, like DNA sequencing, 
gene therapy and cytoplasmic drug delivery in living cells |2|-|5| . Computer simulations have 
been playing a dominant role in most of the expressive literature of polymer translocation, 
which may be, roughly, classified into the topical issues of i) translocation driven by chemi- 



cal potential gradients 
unbiased translocation 
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;10[, ii) translocation driven by external forces 10l-ll3|. and iii) 
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15| . One is usually interested to compute (and, eventually. 



to model) the scaling exponents of the power laws which are found to relate the average 
translocation time r to the polymer size N. 

Unbiased translocation, where the diffusion of a polymer through a membrane pore occurs 
uniquely as the consequence of thermal fluctuations, is by far the most studied case and it 
is also our focus in this paper. Chuang, Kantor and Kardar (CKK) [16| have introduced a 
successful description of unbiased homopolymer translocation from the simple assumption 
that the polymer's evolution does not take it far from its equilibrium states. The essential 
physical picture is that unbiased translocation is ultimately due to the diffusion of the 
polymer center of mass. The translocation time is, thus, assumed to scale precisely in the 
same way as the Rouse relaxation time 17|, so that r ~ A^^"*"^"^, where u is the well-known 
Flory exponent {u ^ 0.588 in three-dimensions). 

Since the scaling behavior r ~ A^i+2^ departs from the one of usual brownian diffusion, 
T oc N'^, it has been eventually suggested [l8|, [l9| that unbiased polymer translocation 
could not be modeled as a markovian process - in other words, memory effects should 
be taken into account as an essential ingredient in any kinematical description of polymer 
translocation. However, it is of crucial importance to note that particular values of the 
translocation exponent are not sufficient per se to rule out the markovian nature of polymer 
translocation. A Markov chain approach can be put forward which actually leads to the 
CKK translocation exponent and to a closed analytical expression for the probability of 
complete polymer translocation that stands in good agreement with results obtained from 
Langevin simulations 20 1. 

Our central aim in this work is to subject the markovian hypothesis of unbiased poly- 



mer translocation to a more stringent test. In essential words, we have considered statistical 
ensembles of three-dimensional polymer translocation realizations, taken from Langevin sim- 
ulations, with sizes considerably larger than the ones presently found in the literature. Our 
statistical data sets have been produced within the Grid Initiatives for e-Science virtual com- 
munities in Europe and Latin America (GISELA), a cloud computing framework supported 
by several academic institutions. 

This paper is organized as follows. In section II, we discuss the Langevin dynamics of 
translocating polymers, which are modeled as bead-spring chains of Lennard- Jones particles 
with the finite-extension nonlinear elastic (FENE) potential. As a starting point, our simu- 
lations are validated from evaluations of the mean translocation time, which is verified to be 



compatible with the CKK scaling prediction 16|. We, then, compute two-point correlators 



of the translocation coordinate and, from them, the memory time of monomer translocation 
events. The polymer diffusion through the membrane pore is noted to be free of anoma- 
lies for times larger than the memory time scale, clearly suggesting that unbiased polymer 
translocation is essentially Markovian. Once we have found, actually, that small monomer 
clusters are correlated during translocation, we take this information into account in section 
III to compare, with good agreement, the empirical probabilities of complete translocation 
for polymers of various sizes to the theoretical probabilities obtained from the Markov chain 



approach of Ref. 20|. In section IV, we summarize our results and point out directions of 
further research. 

II. LANGEVIN SIMULATIONS 

In our Langevin simulations, the excluded- volume and van der Waals interactions between 
beads (monomers and membrane atoms) separated by a distance r are modeled through a 
repulsive Lennard- Jones (LJ) potential with cutoff at length 2^/^cr, where a is the bead 

diameter: 

{Ae[{a/rY^-{a/rf]+e, if r < 2^6^ , 
ULj[r) = i (1) 

I , if r > 2i/V . 

Besides the LJ potential, consecutive monomers are subjet to the Finite-Extension Nonlinear 

Elastic (FENE) potential, 

U,{r) = -\kRlHl-{rlR,f]. (2) 



From the above definition, it is clear that the FENE potential does not allow the distance 
between consecutive monomers to become larger than Rq. 

We have studied polymers with sizes up to 300 monomers, which translocate through 
a pore created by the remotion of a single atom at the center of an 80 x 80 monoatomic 
square lattice membrane. Translocation is dynamically described by the following Langevin 
equations, 

fj f ■ _, rj'T' ■ — * 

m-^ = - E ^n[U,M,) + U,{r,,)] -Q + F,{t) , (3) 

where rjj = |?^j ~ '^j\i ^ is the dissipative constant and Fi(t) is a gaussian stochastic force 
which acts on the monomer with label i, completely defined from the expectation values 

{Fi{t)) = , 

{[h ■ Fi{t)] [n' ■ Fj{t')]) = 2n ■ n'kBnSij5{t - t') . (4) 

Above, n and h' are arbitrary unit vectors, and ks and T are the Boltzmann constant and 
the temperature, respectively. By means of a suitable regularization of the stochastic force, 
we have implemented a fourth-order Runge-Kutta scheme for the numerical simulation of 
the Langevin Equations (|3]). Our simulation parameters are: e = 1.0, a = 1.0 (a is also 
identified with the membrane lattice parameter), ^ = 0.7, k = 7e/o"^, Rq = 2a, ksT = 1.2e. 



The simulation time step is taken to be 3 x lO'^tu, where t^j = ^Jvaa"^ je is the usual 
Lenard- Jones time scale. In the most general case, the initial configuration of the polymer 
has n monomers on the trans-side of the membrane and A^ — n on the cis-side. Translocation 
is allowed to start only after an initial stage of thermal equilibrium is reached for the whole 
polymer. 

As a preliminary validation step, we have checked if the translocation samples produced 
from the Langevin simulations would lead, in fact, to mean translocation times that scale 
with the polymer size as expected on the grounds of the CKK phenomenological theory, i.e., 
r ~ A^" with a = 1 + 2z/. The polymers are initially prepared to be in thermal equilibrium 
with an equal number of monomers on both sides of the membrane. Mean translocation 
times have been obtained from averages taken over ensembles of 70 complete translocation 
processes for each given polymer size. The mean translocation time as a function of the 
polymer size is shown in Fig. [H which in fact indicates a reasonable agreement with the 
CKK scaling exponent. Being confident on the Langevin simulational scheme, we are now 



ready to move to the study of more subtle aspects of polymer translocation. 
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FIG. 1: The mean translocation time as a function of polymer size for unbiased translocation. We 
find the translocation exponent a = 2.17 it 0.06, which is compatible with the CKK prediction, 
a = 1 + 21/, with 1/ = 0.588 (numerical value of the three-dimensional Flory exponent). 



Correlation Effects 

An important issue addressed in studies of polymer translocation refers to the role of 
correlations between monomer translocation events. As we show below, relying upon clear 
numerical evidence, such correlations are short-ranged in time, a fact that paves the way for 
a markovian description of unbiased translocation. 

When a given dynamical system is claimed to evolve in time as a markovian stochastic 
process, it is of course tacitly assumed that the states of the system have been picked up in 
time intervals which are larger than some time scale associated to memory effects. In the 
polymer translocation context, it is convenient to define the dynamical state by the integer 
number s(t), the "translocation coordinate", that represents the number of monomers that 
have crossed the membrane to one of its sides, up to time t (which, by the way, is also treated 
as a discrete variable, since we measure it in units of the Lennard- Jones time scale). The 
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FIG. 2: The dashed hne is the empirical evaluation of G(A) as given by the truncated version of 
Eq. ([5|) with finite M and T. The solid line is the excellent fit provided by exp(— A/9300). 



memory time scale can be defined, in principle, from the decaying profile of the normalized 
two-point correlation function 



G{A) 



lim 



^M v^T 



(5) 



^'^^- E;=iE.=oh(^)]^ 

where Si{t) denotes the i-th sample taken from the ensemble of translocation coordinate 
time series. We have evaluated the right-hand side of ([5]) for polymers composed of A^ = 300 
monomers, in an ensemble of M = 1000 samples, with time bound T = 1.5 x 10^. The initial 
polymer configurations are in thermal equilibrium with 150 monomers on each side of the 
membrane. As it is depicted in Fig. |2l the two-point correlation function ([5]) follows, to very 
good approximation, the simple exponential law G{A) = exp(— cA), with decaying time 
parameter 1/c = 9.3 x 10^. In concrete terms, this result means that after around 9.3 x 10^ 
Langevin iterations, translocation looses memory of the past states. It is interesting, having 
in mind Markov modeling, to find the typical size of correlated monomer clusters that 
translocate within the memory time scale. As we discuss below, this can be achieved from 
an analysis of the diffusive behavior of polymer translocation. 
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The normal diffusive regime of unbiased translocation 

Still considering polymers of size A^ = 300 in three dimensions, the variance of the 
translocation coordinate, that is, ([s(t)]^), is plotted in Fig. [31 Averages are now taken 
over ensembles of 10^ samples. It is clear from that picture that for times larger than the 
memory time scale, {[s{t)Y) is essentially a linear function of time or, in equivalent words, 
monomers diffuse in a normal way through the membrane pore. 
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FIG. 3: Plot of the variance of the translocation coordinate as function of time for N=300. A 
linear regime effectively holds for t ^ 8 x 10^, after an initial period of nonlinear transient behav- 
ior. Statistical averaging is performed in an ensemble of 10^ independent polymer translocation 
realizations. 



It is important to emphasize that our data provides a strong objection to the previously 
predicted anomalous scaling ([s(t)]^) ~ ti+a-^ [16|, which must be viewed now as a misleading 
result derived within the CKK theory of the mean translocation time. Fig. [3] also indicates 



that around the memory time scale t = 9.3 x 10"^ we have \/([s(t)]^) ^ 6.5, which means that 
the correlated monomer clusters are composed of approximately six monomers (we choose 
to round the size of the monomer clusters to six and not to seven, once we observe that the 
crossover to the linear regime in Fig. [3] takes place a little before the memory time scale 



set by the two-point correlation function ([5])). The size of such correlated monomer clusters 
is a crucial ingredient in the Markov chain approach to polymer translocation: if the origi- 
nal polymer is replaced by a "monomer-clustered" polymer (whose size is the original size 
divided by the size of correlated monomer clusters), then the translocation of monomer clus- 
ters, rather than individual monomers, is assumed to generate a truly markovian stochastic 
process. 

III. EVIDENCE OF MARKOVIAN BEHAVIOR 

The fundamental hypothesis of the Markov chain approach to unbiased translocation is 
that monomers translocate in an uncorrelated way with probabilities Pn and g„ for cis— J-trans 
and trans— j-cis transitions, respectively (recall that n, as defined in Sec. II, is the number 



of monomers on the trans-side of the membrane). Following Ref. 20[, one puts forward the 
transition probabilities 



Pn 



{N - n)'5+2'^-l 

c 

,(5+2i/-l 



In = ;^5^^ , (6) 



where < c < 1 is an arbitrary constant and 5 is a scaling exponent associated to finite-size 
corrections to the CKK scaling relations (it would follow, for instance, that r ~ ]\[^+'^'^ for 
fixed 6 and large enough A^). It is also possible, from the definitions ([6]) and using exact 



results for general Markov chains 21[, to find the probability P{N, n) of complete trans— !>cis 
polymer translocation, 

P(Ar, n) = ^^=\^^^='P. . (7) 

i + E^^^nkS 

We have tested the analytical prediction given by Eq. ([7]) for polymers of various sizes. 
As it is shown in Fig. |U the comparison between the empirical and analytical probabilities 
is very satisfactory. The empirical probabilities have been evaluated from ensembles of 70 
complete polymer translocation realizations. To understand Fig. HJ note that the parameter 
A^ to be substituted in ([6]) is not the original polymer size anymore. In ([6]), A^ is now taken 
as the effective size of the monomer-clustered polymer (in our particular case, as suggested 
by the results of Sec. II, it is just the original polymer size divided by 6). Therefore, for 
a given value n of the initial number of trans-monomers, a solid circle is plotted in Fig. HJ 
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n: Initial Number of Trans-Monomers 

FIG. 4: The empirical (open circles) and analytical (solid circles) probabilities (both denoted here 
by P{N,n)) of cis — > trans complete translocations are compared. The variable n stands for the 
initial number of monomers in the trans-side of the membrane. The analytical probabilities are com- 
puted from dl]) and (HD with 5 = 0.88 and u = 0.588 for N = 60, 70, 90, 100. The theoretical error 
bars follow from elementary statistical considerations and are given by ^yP{N,n)(l — P{N,n))/70. 



with coordinates (n, P{N/6, n/6)) where P{N/6, n/6) is the probability of complete polymer 
translocation evaluated within the Markov chain approach for a polymer which contains N/6 
monomers (A^/6 is conventionally rounded, if necessary to the smallest integer greater than 

A^/6) . 



IV. CONCLUSIONS 

We have provided consistent statistical data which essentially settles down the issue on 
whether unbiased polymer translocation is markovian or not - it turns out that unbiased 
polymer translocation can be very confidently described as a markovian stochastic process. 
This conclusion is supported from three clear pieces of evidence: (i) the two-point correlation 
function ([5]) indicates that individual monomer translocation events are correlated whithin 
a finite memory time scale which is much smaller than the mean complete translocation 
time; (ii) the variance of the translocation coordinate depends linearly on time for times 
which are larger than the memory time scale; and (iii) empirical probabilities of complete 
translocation finely match the analytical ones predicted from the Markov chain approach of 



Ref. 20| (it is worth of mentioning that (ii) can be also derived within the same formalism). 



An interesting point, motivated by our results, is whether the original approach of poly- 



mer translocation addressed by Muthukumar 22| , which is also markovian (but failed to give 



the correct expression for the mean translocation time), can be somehow improved taking 
into account the present findings. In concrete terms, it is not difficult to define a specific 
free-energy profile for the translocating polymer that would lead to the individual monomer 
translocation probabilities (|6]). A challeging problem which we are currently investigating 
is whether such effective free-energy profile is just an artifact that can be used to reproduce 
the stochastic evolution of the translocating polymer, or is actually the physical thermody- 
namical potential derived from standard equilibrium statistical mechanics considerations. 

As a final remark, we call attention to recent works on driven translocation, where memory 
effects are found to be relevant, due to the excitation of collective modes along the polymer 
chain 



11 



13| . However, as a problem deserved for further research, we note that if the 
forces which pull the polymer are not strong enough, it is likely that translocation can still 
be described with the help of Markov chain equations, where the transition probabilities ([6]) 
are just replaced by alternative ones. 

This work has been partially supported by CNPq and FAPERJ. 
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